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Based on the discrete fractional random transform (DFRNT), we present the dis- 
crete fractional random cosine and sine transforms (DFRNCT and DFRNST). We 
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1 Introduction 



Recently we proposed a discrete fractional random transform (DFRNT) [1] and 
have demonstrated its application to image encryption and decryption. The 
DFRNT is a kind of discrete transform with fractional order originated from 
the fractional fourier transform (FrFT) [2] and especially the discrete frac- 
tional Fourier transform (DFrFT)[3], and thus has the most excellent mathe- 
matical properties as FrFT and DFrFT have. Meanwhile, however, the result 
of transform itself can be inherently random. 

Because the close relationship of the DFRNT and the DFrFT(FrFT), we hope 
that the DFRNT can serve as a general mathematical and numerical tool in 
the field of digital signal processing. In the first step towards this end, we 
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extend the DFRNT to discrete fractional random cosine and sine transforms 
(DFRNCT and DFRNST). We know that the Cosine and Sine transforms and 
their discrete versions are useful tools in signal and image processing, such 
as signal coding[4], watermarking [5] and restoration of de-focused images[6]. 
In the Ref. [7] Pei and Yeh extended the Cosine transform to the discrete 
fractional cosine transform (DFrCT) and the discrete fractional sine trans- 
form (DFrST). Both of them possess well the angle additivity property of the 
DFrFT. Moreover, the DFrCT and DFrST are used in the digital computation 
of FrFT for the reducing computational load of the DFrFT. 

The motivation of this paper is to propose the definitions of the DFRNCT 
and DFRNST directly based on the DFRNT and investigate their mathe- 
matical properties. We find that the DFRNCT and DFRNST, generated by 
the method of Pei and Yeh [7], are very close to the original DFRNT. We 
can demonstrate that they are only two subsets of DFRNT. From DFRNCT 
and DFRNST we can also regenerate another DFRNT which will have nice 
symmetric properties for even and odd signals. Numerical simulations have 
demonstrated such properties. 

In Section 2, we briefly introduce the definition of DFRNT and give its 
mathematical properties. And then we give the definitions of DFRNCT and 
DFRNST based on the DFRNT. Their mathematical properties and the re- 
lationships with DFRNT are also given. In Section 3 we give some numerical 
results. Conclusions are given in Section 4. 



2 The DFRNT, DFRNCT and DFRNST 

2. 1 The DFRNT and its properties 

The DFRNT's of a 1-D and a 2-D signals, denoted by x and y respectively, 
can be written as matrix multiplications as follows [1] 



where lZ a is the kernel transform matrix of the DFRNT and can be expressed 
as 



X a = TZ a x, 



(1) 



Y a = K a y{K a )\ 



(2) 



n a = vD a v\ 



(3) 
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In the kernel matrix, D a is a diagonal matrix generated by a set of values 
{exp(— 2iima/M) : n = 0,1,2, ...N — 1}, which are considered to be the 
eigenvalues of the DFRNT. Where a indicates the fractional order of the 
DFRNT. M is a positive number, usually is an integer which has a meaning 
of periodicity with respective to the fractional order a in eigenvalues. D a is 
written as follows 



D a = diag 



2ma Airca. . 2i(N — l)ira 
1, «p(— i^), exp(-— ), exp( ) 



•(4) 



The randomness of the transform comes from the matrices V and V 1 , where 
[. . .]* indicates the transpose matrix of the matrix [...]. The matrix V is gen- 
erated by iV orthogonal vectors {vi,v 2 , ■■■,vn} as 

V = [vi, v 2 , v N }. (5) 



Where v n (n = 1,2,...N) are column vectors which are the eigenvectors of a 
symmetric random matrix Q with Q mn = Q nm - The matrix Q can be obtained 
by an N x iV real random matrix P with a relation as 

Q=(P + P t )/2. (6) 



The matrix V satisfy the relation of VV* = /. 

Even though the DFRNT results in a random vector, the transform itself still 
has several excellent mathematical properties as follows. 

[1] Linearity. ll a {ax x + bx 2 ) = dR a x x + bTZ a x 2 . 

[2] Unitarity. TZ' = (K a )*. 

[3] Index additivity. ll a Vf = Vfll a = K a+ P . 

[4] Multiplicity. U a+M = Tl a . 

[5] Parseval energy conservation theorem. J2 \X a (K)\ 2 = J2 \x(k)\ 2 . 

K k 

2.2 DFRNCT, DFRNST and their properties 



The definitions of DFRNCT and DFRNST can be directly given from the 
DFRNT. Similar to the definitions of DFrCT and DFrST, of which the eigen- 
vectors can obtained from the eigenvectors of the DFrFT[7], the eigenvectors 
of DFRNCT and DFRNST are the same with that of the DFRNT if the same 
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symmetric random matrix Q are used. That means when we construct the 
DFRNCT and DFRNST, we use the same method to generate the matrix 
V = [vi, v 2 , . . • , vjy}. However, in the following discussion, we denote the the 
eigenvectors of the DFRNCT and the DFRNST by matrices 

V c = [Ci, C 2 , ... , Cat], V s = [si, s 2 , . . . , s N }. (7) 



Bear in mind that c n = s n = v n for the same matrix Q. 

Similarly the kernel matrices of the DFRNCT and the DFRNST are defined 
as follows 

R a c =V c D a c (V c )\ R" = V s D%(V s y. (8) 



However the two eigenvalue diagonal matrices for DFRNCT and DFRNST, 
i.e. and Df are, respectively, chosen as follows 



D a c = diag 
Dg = diag 



Aina, , 2(2N -2)ma. 
1, ex P( — ex P( V} ) 



exp( 



M 
2ina 



M 



M 



Qiira 2(2N — l)ina 
h ex P( Jfh •••» ex P( u ) 



M 



(9) 
.(10) 



The eigenvalue matrices of the DFRNCT and DFRNST are also diagonal ma- 
trices with the diagonal elements chosen within the set of {exp(— 2inna/M) : 
n — 0, 1, 2, ...2N — 1}. The elements of the DFRNCT are formed by the values 
with n = 2j,j = 0, 1, 2, . . . while the elements of the DFRNST are given by 
the values with n — 2j + 1. The parameter a indicates the fractional order of 
the transform and M relate with the periodicity. 

The DFRNCT and DFRNST of a 1-D signal x and a 2-D image y then can 
be expressed as 

x: = r?x, x: = r: x , (n) 

Y? = RZy[RZ\*, Y s a = R:y[R:} t . (12) 



When a — 0, the kernels of the DFRNCT and DFRNST are identity matrices. 
The DFRNCT is an identity matrix for a = M/2, however the DFRNST is a 
negative identity matrix (Rfl 2 = —I) for a = M/2. 

From the Eq. 9 and Eq. 10 one can see that the DFRNCT and DFRNST 
defined here are not generated by merely taking the real and imaginary parts 
of the DFRNT kernel matrix. Therefore, the construction of DFRNCT and 
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DFRNST are different from the FrCT and FrST in the Ref . [8] . If we define the 
DFRNCT and DFRNST that way, they would lose the property of additivity, 
and they would be a discrete transform in real number domain. And more 
significantly in such cases the corresponding inverse transforms do not exist. 

On the definitions of DFRNCT and DFRNST, we adopt the method of Pei 
and Yeh[7]. Because we use orthogonal eigenvectors to define the DFRNCT 
and DFRNST, the mathematical properties of DFRNCT and DFRNST are 
thus similar to DFRNT. 

• Unitarity. This property is directly from the diagonal matrix (and Df), 
that is 

[R a r = [RT 1 = R -a^ (13) 
[R * r = [RT 1 = R -a (M) 

• Index additivity. This property always valid for DFRNCT and DFRNST 

because 

R2RS = RgRZ = R^. (15) 

iffitf = RgRg = R^. (16) 

Moreover, according to this property, the matrices of inverse transform of 
DFRNCT and DFRNST with fractional order a are R~ a and R; a , respec- 
tively. 

• Periodicity. The DFRNST is periodic with M, however the DFRNCT is 
periodic with M/2, i.e. 

R a+M/2 = R a^ R a+M = R a ^ 

• Energy conservation theorem. 

EI^M 2 = EI^M 2 = EI^)I 2 - ( 18 ) 



2.3 The relationships between DFRNCT, DFRNST and DFRNT 



The general relationships between DFRNT, DFRNCT and DFRNST can be 
explored from the relationships of their eigenvalue matrices, which are given 
by the equations Eq. 4, Eq. 9 and Eq. 10. All the transforms take the same 
eigenvector matrices, the only differences are their eigenvalues. From the def- 
initions, one can find immediately that 

R a c = TZ 2a , R a s = exp(-2ia7r/M)n 2a . (19) 
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Therefore, the DFRNCT and DFRNST can be regarded as DFRNT's with 
changes of scale in the fractional orders. In this instance, the DFRNCT and 
DFRNST form two subsets within the whole domain of DFRNT. 



In above discussion we assume that all the transform kernels have the same 
N x N dimension. We can also regenerate a DFRNT with all the eigenvectors 
and eigenvalues of DFRNCT and DFRNST. Such reconstructed DFRNT will 
have a 2N x 2N or (2iV + 1) x (2iV + 1) transform kernel matrices which 
can process the data with 2N- and (2N + l)-points long. In order to indicate 
the difference between the reconstructed DFRNT and the original definition, 
we refer to it as ReDFRNT. We assign the following eigenvalue matrices to a 
ReDFRNT 



D 



2N 



diag 



1, exp( 



2ma, 
M ' 



cxp( 



2{2N - l)ma. 

M ' 



(20) 



and 



D 



27V+1 



diag 



1, exp( 



2ma 
M 



), exp( 



ANina 
M 



(21) 



For the 2iV-point ReDFRNT, the eigenvectors is given by the eigenvectors of 
DFRNCT and DFRNST as the following matrix 



Von — 



V2 



Cl, Si, C2, S2, 



, Cat, Sat 



(22) 



Where c z n (or s^, n 



1,2, ...,iV) denotes the flipping eigenvector c n (or s n ) 
in up-down direction. The eigenvectors V2N are orthogonally vectors, i.e. 
V2nV2 N = I. For a (2iV+l)-point ReDFRNT, however, the eigenvectors should 
be expressed as 



V2N+I — —F= 



1 

V2 



Cl, Si, c 2 , s 2 , 
0, 0, 0, 0, 



, C N , S N , V 

, 0, 0, V2 

5 c Ni ~ s Ni v 



(23) 



so that V^at+i can also be an orthogonally eigenvectors, where vo is zero vector. 
Then the kernel matrix of the ReDFRNT can be written as 



(24) 



where m denotes 2N or 2N + 1, respectively. 
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Obviously the ReDFRNT is different from the DFRNT originally defined. Here 
the ReDFRNT, similar to the DFRNCT and DFRNST, can be regarded as a 
subset of DFRNT. Such relationships resemble the self-affine characteristics 
of fractals. 

The eigenvector matrices of the ReDFRNT can also be written as the format 
of column vectors, for instance V 2 n — [vi v 2 ... v 2 n] , where 









Sn 


V2n-1 = 




, v 2n = 






c z 




— S 2 



The eigenvectors v 2n -i and v 2n are even and odd eigenvectors which have even 
and odd symmetry, respectively. Same results are also valid for the matrix 
V2AT+1- Because the eigenvectors of the ReDFRNT have well symmetry, the 
outputs can be calculated with DFRNCT and DFRNST for even and odd 
signals, respectively. The amplitudes of the ReDFRNT have even symmetry 
for both of even and odd signals, and its phases have even symmetry for the 
even signal. However, for the odd signal the symmetric relationship of phases 
is not obvious because 

<f>[n) = <p(2N + 1 - n) ± tt. (26) 



Where 4>(ri) is the phase of the ReDFRNT, n can be taken the value from 1 
to 2N. We then define a special phase for the ReDFRNT of an odd signal as 

0'(n) = arctan{tan[0(n)]}, <j>'(n) G [-tt/2, tt/2], (27) 



Such defined special phase 4>'{n) will be symmetric for an odd signal, as de- 
picted in Fig. 3. 

The above proposition can be proved from the symmetry of the eigenvectors 
of the ReDFRNT. Numerical calculations can also demonstrate this proposi- 
tion (see Fig. 1 to Fig. 3). The symmetrical properties will be useful in the 
practical computation of the reconstructed DFRNT. The eigenvectors of the 
DFRNCT and DFRNST can be chosen as uniformly distributed values, thus 
the computational load can be reduced to one half. 

We know that in the continuous function f(x) can be decomposed into 

an even and an odd functions(/ e (a;) and f (x)) as 

Ux) = [/(«) + /(-*)] , Ux) = [/(*)-/(-*)] . (28) 
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Similarly a discrete 2iV-point signal s(n) can also be decomposed into an even 
and an odd signals(s e and s ) as 



(29) 
(30) 



Because the amplitude and the phase of a ReDFRNT have well symmetry 
for even and odd signals, we only need to compute the DFRNCT and the 
DFRNST of half of signal(s' e (n) = s e (n) and s' (n) = s (n), n = 1, 2, N) as 



a _ d«„' c _ pa r 



(31) 



And the output S'(n) of the ReDFRNT can be computed with S ec and S os as 

S ec (n)+S os (n) if l<n<iV, 

5(n) = { (32) 

5 ec (2A^ + 1 - n) - S os (2N + 1 - n) ii N < n < 2N, 



For (2iV + l)-point signal s(n), the output of ReDFRNT can been expressed 
this equation as follows 



S(n) = 



S ec (n) + S os (n) 
s(N + l)exp(- ANi7Ta * 



M 



if 1 < n < N, 
if n = N+ 1, 



(33) 



,S ec (2A^ + 2 - n) - S os {2N + 2 - n) if TV + 1< n < 2N + 1. 



Thereby, we can use an iV x N kernel matrix to count the ReDFRNT of an 
arbitrary 2iV-point(or (2N + l)-point) discrete signal. 



3 Numerical results of simulation 



We know that the ReDFRNT, restructured by the DFRNCT and DFRNST, 
provides symmetric distribution for both even and odd signals. Such property 
can be easily demonstrated by the following numerical calculations. 

In the computation we choose the following two simple functions xi{n) and 
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X2(n) with size of 128 points, 



xi(n) 



1 



if 49 < n < 80, 



(34) 



< 







otherwise, 



and 



1 



if 49 < n < 64, 



x 2 {n) 



< 



-1 



if 65 < n < 80, 



(35) 







otherwise. 



Where x\{n) and x 2 in) are even and odd symmetric signals, respectively, 
with 1 < n < 128. The random matrices P c and P s are generated with 
MATLAB, and the fractional order a and M are chosen as a = 0.6 and 
M — 1, respectively. In order to compare the numerical results of DFRNCT 
and DFRNST with ReDFRNT, we only take one half of the above signal X\ 
and X2 for DFRNCT and DFRNST, i.e. the signals x[ = Xi{n) and x' 2 = X2(n), 
n = 1,2,. .,64, respectively. The corresponding eigenvalues are calculated ac- 
cording to Eq. 9 and Eq. 10 with N = 64. The matrices of the eigenvectors are 
computed from the random matrix Q (Eq. 6), indicated by Eq. 5 and Eq. 7. 
The ReDFRNT's of the signals X\ and x 2 are calculated from the eigenvalues 
given by Eq. 20 and Eq. 21, with the corresponding eigenvector matrices given 
by Eq. 22 and Eq. 23, respectively. 

The results of one dimension DFRNCT and the DFRNST are given in Fig. 1 
and Fig. 2. The bold lines denote the amplitude of the DFRNCT and DFRNST 
for the half signals x[ and x' 2 , which coincide well with the ReDFRNT's am- 
plitudes of X\ and x 2 , respectively, when n < 64. The symmetric properties 
of ReDFRNT for an even signal in both amplitudes and phases have been 
revealed clearly in the figures. The phase of the ReDFRNT for signal x 2 is 
not symmetric, however the special phase 0'(n) defined by Eq. 27 can be to- 
tally symmetric for the ReDFRNT of the odd signal x 2 . Such result is shown 
in Fig. 3. Where we can find both amplitude and phase 4>'{n) are symmetric 
with respect to the central line x = 64.5. 

For the case of two dimensional transforms, three binary images with 128 x 128 
pixels Ji, J 2 and J 3 , shown in Fig. 4(a)-(c) respectively, are used. The image 
only contain simple rectangular patterns which is equivalent to some rectan- 
gular functions of rect(a; — xo, y — yo)- Fig. 4(d)-(f) illustrate the corresponding 
results of ReDFRNT. Here only amplitudes are displayed. The same param- 
eters are chosen as the case of one dimensional transforms. The results show 
that the ReDFRNT keeps the symmetry of input images. 
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4 Conclusion 



In this paper, we proposed the discrete fractional random cosine and sine 
transforms (DFRNCT and DFRNST) based on the discrete fractional ran- 
dom transform (DFRNT). These two random transforms are generated by 
rearranging the distributions of eigenvalue matrices meanwhile keeping the 
eigenvector matrix unchanged. Such defined DFRNCT and DFRNST have 
excellent mathematical properties inherited from DFRNT. We have demon- 
strated that the DFRNCT and DFRNST are nothing but scaled DFRNT and 
form two subsets of DFRNT. From these two transforms we can also recon- 
struct another kind of DFRNT (ReDFRNT) which combines DFRNCT and 
DFRNST together. The ReDFRNT has a special feature that it has symmetric 
distributions for even and odd signals in both amplitudes and phases. 

The DFRNT, DFRNCT, DFRNST and ReDFRNT are discrete fractional or- 
der transforms with intrinsic randomness. We have demonstrated that the 
DFRNT can be applied in information security such as image encryption and 
decryption. Further applications of these random transforms in image pro- 
cessing, pattern recognition, artificial intelligence and so on are left as open 
questions for the community of information science. 
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List of figure captions 



Figure 1. The results of DFRNCT and ReDFRNT of the 1-D signal x u with 
a = 0.6 and M — 1. Because the function x\ is even, so that the ReDFRNT 
are symmetrical and coincides with the DFRNCT. 



Figure 2. The results of DFRNST and ReDFRNT for the odd signal %2 with 
a = 0.6 and M — 1. The amplitude of ReDFRNT are symmetric and coincides 
with the result of DFRNST. 



Figure 3. The same calculation results of DFRNST and ReDFRNT for the 
signal x 2 with special phase <t>'{n). Where we also take a = 0.6 and M = 1. 
The special phase then is symmetric distributed for the ReDFRNT. 



Figure 4. The results of ReDFRNT for 2-D functions: (a), (b) and (c) dis- 
play three input images ii, I 2 and I3, respectively. From (d) to (f) are the 
corresponding ReDFRNT's, with a = 0.6 and M — 1, for images Ji, I2 and 
I3, respectively. From the results we can find that the ReDFRNT keep the 
symmetries in the input images. 
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Fig. 1. The results of DFRNCT and ReDFRNT of the 1-D signal x 1 , with a = 0.6 
and M = 1. Because the function x\ is even, so that the ReDFRNT are symmetrical 
and coincides with the DFRNCT. 
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Fig. 2. The results of DFRNST and ReDFRNT for the odd signal x 2 with a = 0.6 
and M = 1. The amplitude of ReDFRNT are symmetric and coincides with the 
result of DFRNST. 
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Fig. 3. The same calculation results of DFRNST and ReDFRNT for the signal X2 
with special phase (f>'(n). Where we also take a = 0.6 and M = 1. The special phase 
then is symmetric distributed for the ReDFRNT. 
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(d) (e) (f) 

Fig. 4. The results of ReDFRNT for 2-D functions: (a), (b) and (c) display three 
input images ii, I2 and ^3, respectively. From (d) to (f) are the corresponding 
RcDFRNT's, with a = 0.6 and M = 1, for images I\, I2 and ^3, respectively. From 
the results we can find that the ReDFRNT keep the symmetries in the input images. 
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